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ABSTRACT j 

The Airborne Laser Ranging System is a proposed multibeam short pulse laser ranging system i 

on board an aircraft. It simultaneously measures the distances between the aircraft and six laser : 

retroreflcctors (targets) deployed on the Earth’s surface (Figure 1). The system can interrogate 

over 100 targets distributed overan area of 2.5 X Iff* sq. kilometers in a matterof hours. Poten- f 

tially, a total of 1.3 million individual range measurements can be made in a six hour flight. The 
precision of these range measurements is approximately ± 1 cm ( 1). These measurements arc tlicn 

used in a procedure which is basically an extension of trilateration techniques to derive the intersite 
vector between the laser ground targets. By repeating the estimation of the intersite vector, strain 
and strain rate errors can be estimated. These quantities arc essential for crustal dynamic studies 
wliicli include determination and monitoring of regional strain in the \icinity of active fault zones, 
land subsidence, and edifice building preceding volcanic eruptions. 








Till: AIRBORNH LASliR RANGING SYSTEM, 

ITS CAPABILITIBS AND AITLICA TIONS 

1.0 INTRODUCTION 

Recent experience with laser and Very Long Baseline Interreromeler (VLBI) measurements in 
Southern California reveal that large scale crustal motions can occur in time scales of a few weeks 
and months. To adequately monitor such motion, techniques are required to map the position of 
grid points over a region in a few days and at repeat frequencies of a few weeks. In addition, 
maps of crustal dcforniation rates encompassing 20 to 40 locations aa* aniuired within several 
days of the observations. The monitoring of relative motions in the Barth’s upper crust in tens 
to sub-centimeter rale per year can be accomplished by a ptilsed laser ranging system carried on- 
board an aircraft, making rapid range measurements to passive reHeclors distributed on the ground. 
By developing and interpreting this system's ability to delect motions of the Barth’s upper crust, 
a model of the strain accumulation compatible with observations of crustal motion and tectonics 
of a region within the experimental data collection area can be derived. B'urlhennore, an Air- 
borne Laser Ranging System (ALRS) can sur\ey an area in a very short period of time (hrs.) and 
resurvey the areas as required. 

The basic philosophy of the Airborne Laser Ranging System is to invert the usual laser rang- 
ing configuratiofi by placing the ranging and pointing hardware in an aircraft such as NASA’s 
NB3A Lockheed Orion Research Aircraft, and replacing the expensive laser ground stations by 
low cost (< SI 000) passive retro re Hectors. The system is necessarily multibeam since the loca- 
tion of the aircraft is not known with cm precision at each point where a set of range measure- 
ments are made. Thus, a minimum of four simultaneous range measurements are required, i.c., 
three to resolve the new coordinates of the aircraft and one to acquire information on the rela- 
tive locations of the ground targets. The ALRS system will he capable of ranging simultaneously 
to six retrorelleclors. At a l.iser repetition rale of lOpps, a potential 1.3 million individual range 
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an b. ,n,.c anO an ,«a as ,a«n as 60.000 s,. km can b. sumcycO dn„n,onc s,x 
bon, night. Tb. latte, cove, age arpbes to a high altittnlc mseateh ai,e,aft such as an RO-57 «, - 

altitudes. 

compote, simulations have shown that in the pmsenee of measumment noise and bias ensupled 
with tmpospbetle tefmetion effeets. an ai,e,aft oiwmting at a mote modest maximum altitude of ^ 
6km. can detennine intentite distances to a ptecision of 0.4e,n at 5 km and 1.4em at 30km base me 
distaneew. The e„o, gmwth mte pe, unit baseUne vanes invemely with ahetaft altitude. Futthc 
mom the data mduetion teehnhtue simultaneously msolves the ai, craft position to the cm level at 
each point in the night path whem a Ime, pulse is ttansmitted. The AbRS system is expeetexl to be 
a powetful new ,eseateh too, fo, monitodug mgiona, costal motion, land mauagement appheattons. 
and gene, a, suceying beeaum- „ will pnavlde a -snamhof of the tatget positions ove, an extended 

area with hipli spatial resolution. 

2.0 SYSTl-M DF.SCRirriON 

2.1 Laser Ranging Subsystem 

Flgum : is a block diagtam of the ALRS. The system compute, enables the firing of a sub- 
nanosecond lase, „an.mi„e, af a nominal ,ate of lOpPs. The ttansmitte, is a modeloekvnl. ITM 
Q,w„ehed Nd:VAC, lase, oseillato, followed by a double-pass Nd:VAO lase, amplifte, and a 
Kt,>P netiueney double,. On each fulng, the riansmitte, gencates a single ISOpsee (FWllM) 
pulse containing seveal milliioules of enemy at the 0.53: miemmete, green wavelength. A beam- 
,p„„e, ,eneets a ve,y small action ,< 1. 1 «f .he outgoing enemy into a setic of six beamspht- 
tes Which divide and dhect the low-level enegy into each of six teeeive, channels. The temain- 
ing enegy is divided into approximately six e.,ua, part, by a second set of bea.u.phtters wbreh 
directs the energy to six independently controlled pointing systems. The six outgoing PUlsc Pass 
ttuough the atmosphere to ate si, ground target retroreneetors. The reOeeted energy 
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I’igurc 2. Block Diagram of the Airborne Laser Ranging System 
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from each target travels back through the atmosphere to the ALRS and is imaged onto the corres- 
ponding high speed photomultiplier tube (PMT). Thus, a pair of start and stop pulses, indicating 
the times at which a given pulse leaves and returns to the instrument, are recorded by each of 
the six receiver channels. The use of common start/stop receiver components eliminates a poten- 
tial source of time-dependent range bias whicli might be introduced by changes in the thermal en- 
vironment, voltage condition, or other nonstationary processes during flight. 

Tlic output of each PMT is split inside the signal processor. One port is input to a low time 
walk constant fraction discriminator. The discriminator provides a NIM logic pulse to both start, 
and later stop, one channel of a multichannel time interval unit. The latter device measures the 
time-of-flighl for each of the six laser pulses. The six channels in the TIU share a common clock 
input. The range to each target is calculated from time-of-flight with suitable corrections being 
made for instrument biases, pulse amplitude effects, and atmospheric refraction delays. A charge 
digitizer at the second port of the PMT output records the energies of the outgoing (START) and 
incoming (STOP) pulses thereby enabling the system computer to correct for timing biases due to 
pulse amplitude (dynamic range) effects. The measurement data is transferred to the system 
computer for storage and use by the navigational subsystem as described in Section 2.3. Nominal 
atmospheric refraction corrections are made inflight for use by the navigation computer. More 
exact corrections are made during the post-flight data reduction phase. 

Instrument related single-shot range uncertainties are about 5mm RMS (one sigma) with an 
average received signal level of 200 photoelectrons. This is the signal level calculated for a worst- 
case link which assumes 1 mJ of laser energy per channel, a 0.5° beam divergence, and a min- 
imum aircraft elevation angle of 20° as viewed by a ground target with a modest cross-section of 
lO^m^. The aircraft is assumed to be at its maximum altitude of 6km. 
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The laser transmitter, beam splitting optics, receiver optics, anil photomultiplier tubes arc 
mounteil on an optical baseplate which is isolated vibrationally from the aircraft fuselage. Six 
azimuth-elevatio!! pointing mounts with 5 cm n.'tvive apertures arc rigidly attached to the bottom 
of the optical bed. I-ach pointing system c-onsists of a four mirror coclostat mounted on an 
azimuthally rotating stage as in Figure 3. The laser beam enters the pointing system through a 
hole in the optical b.iseplate and the rotation stage. The final mirror rotates about an axis paral- 
lel to the optical bench to point to a given elevation angle. This particular configuration was 
chosen because it can be placed ver>' close to the aircraft window to provide near hemispherical 


FIXED 



Figure 3. Optical Tointing System Concept 
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viewing capability. A secoiul advantage is that all light beams reflect at a 45® incidence angle 
independent of the AZ/EL pointing angles. This allows the use of high efficiency dielectric coatings 
on the mirrois Ter maximum receiver sensitirity and potential upgrade of the system to two colors 
for making direct measurements of the atmospheric refraction cora'clion. (I) 

Each pointing mount is eiiuipped with two servo systems (azimuth and elevation). In the 
present concept, these systems are digital controllers that drive precision stepping motors and con* 
tain optical encoders to measure angular position. Tlie commands sent to the controller arc in 
the form of angular position, velocity, and acceleration referred to the optical bench. The con- 
troller performs tire usual loop closure tasks to control the individual servos as directed and to 
relay information back to the command interface. Command angles and predicted range gates 
are generated by the system computer using aircraft navigation solutions provided by the naviga- 
tion and attitude determination subsystem. 

Witli laser beam divergences on the order of several milliradians, absolute pointing accuracies 
at the mrad level are adequate. Tliis perfonnance is a factor 10 to 20 times less stringent than 
typically required for ground-based satellite laser ranging systems. 

2.2 Target Heployment 

A number of geodetic monuments will be erected in the region of interest. The ‘•grapefruit- 
sized” targets, studded with about six optical cube comers, can be mounted directly on the 
monuments prior to the suney flight. These can be cither permanent installations or the targets 
can be removed and reused at other locations. In order to tie the target grid accurately to a 
reference cv>orvlinate s>^tem targets would be placed at three or mon: Hducial points. To sim- 
plify the target aaphsilion sequence, it is desirable that the position of each target be known 
apriori to approximately 50 meters but this is not a hard requirement. The apriori knowledge 
of target position does not impact the aposleriori grid absolution adiieved by the ALRS. In 
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general, apriori location can be read from a siincyors map or obtained using radio navigation 
a-ccivcR such as Loran C, tlie TRANSIT Satellite System, or the futum Global l\)sitioning 
System (GPS). 


2.3 Target Acquisition and Tracking 

The a priori target positions discussed in the previous section will be stored in the ALRS sys- 
tem computer mcmor>’ as latitude, longitude, and altitude above sea level. Successful acquisition 
and tracking of the targets during flight requires an adequate knowledge of the aircraft position, 
velocity, and attitude. Tlie positional error in a modem inertial navigation system (INS) typically 
grows at a rate of 0.4 to 4.0km/hour. In such a system, vehicle accelerations and attitude are 
measured by an inertial measurement unit (IMU) consisting of three orthogonal accelerometers 
and a gyro triad. The navigational computer (NC) performs coordinate transformations and inte- 
grates the equations of motion to provide estimates of velocity and position relative to a set of 
initial a>iulitions. The errors in these estimates have many contributing sources including sensor 
calibration limitations and inaccuracies, computational errors, and sensor-error-propagation effects, 

Tlie performance of an INS can be improved significantly by incorporating additional inde- 
pendent sensors which periodically check and update the navigation solutions. Radio navigation 
aids (such as OMIiGA, LORAN C\ and the future GPS system) are particularly well suited to the 
role of auxiliary sensor because of their global or near global coverage. In the ALRS. latitude 
and longitude information from a LORAN C receiver is utili/ed to update and stabiii/e the cor- 
responding INS solutions via a Kalman filler algorithm in the NC (2). Similarly, barometric altim- 
eter measurements of altitude are processed in the NC to stabilize the equations for vertical velo- 
city and altitude. Tlie Kalman filter also updates attitude information and pnnides best estimates 
of sensor errors such as misalignments, acccicromctcr biases and scale factors, gyro drift rales, etc. 
The INS velocity solution .issisls the LOR.AN C receiver in the acquisition of Doppicr-shifted 
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signals from the LORAN ground network. In the event of certain types of failure, the proposed 
navigational system can be restarted and calibrated in fliglit. 

Fhglit tests of a LORAN-aided INS (3) perfonned for the U.S. Air Force have demonstrated 
a 60 meter absolute position accuracy (one sigma) and a 15 meter RMS random noise error (one 
sigma) and angular accuracies of a few tenths of a niilliradian. LORAN C was chosen for the 
ALRS because it is significantly more accurate than OMFGA (±925 meters) and, unlike GPS. is ‘ 
already available in most regions of the world. The Global Positioning System will eventually 
provide more accurate dynamic fixes in all three coordinates (5 mctcis) and global coverage, but 
the system is unlikely to become fully operational before 1988. 

Tlie navigational data is combined with the stored a priori target positions to compute the 
estimated range gates and the pointing mount command angles. As the targets come within range, 
the laser is activated and the presence or absence of range returns is noted. If no returns are de- 
tected. a search pattern is executed until range data is acquired. Selected range data is then 
passed to the navigation Kalman filter to update the estimate of aircraft position. The ALRS 
then shifts from the acquisition to the tracking mode. 

In the tracking mode, triangulation on the highly accurate laser returns results in an aircraft 
position estimate which is virtually error free (better than a meter standard deviation) so that, 
except for aircraft attitude estimation errors contributed by the gyros, the computed command 
angles and range gates are essentially correct. The instrument remains in the tracking morle as 
long as sufficient range data is available to maintain an accurate estimate of aircraft position. In 
making the transition to a new set of six targets, the system triangulates on laser returns which 
arc common to the new and the previous set. In this way. the aircraft position is known with 
better than meter accuracy during the transition. The design of the system computer and a 
description of the inflight operational software is given elsewhere. (4) 


8 









.^ ' ' . J * :- ^' ■J^;^T^A ■ V' » **/ * ' - *" ”-^y. 


3.0 MATIIOMATICAL MODELS AND ALGORITHMS 
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The mathematical model associated with the data generated by the ALRS system will be 
given in this section. This model will lead to the definition of an estimation algorithm for the 
coordinates of the retroreflector positions. 


3.1 Range Data Model 

Ideally, the ranging system measures the distance (range) between the airborne Ixser and 
several retroreflcctors deployed on the Earth’s surface. 


2(3X 1)(0 iif'J ^(3Xt)(0 be the 3-dimensional vectors describing the position and velo- 
city of the aircraft at time t. 


^3Xt)t*) " 

Z|(t) 

Zs(t) 

and Z(jx,)(t) = 

Z|(t) 

ij(t) 


Z 3 O) 

(3X1) 

Z 3 (t) 


The vector components are expressed in some convenient, coordinate system. 


( 10 ) 


Also let, 

- 

be the 3-diniensional vector describing the position of the retroreflector in the same coordi- 
nate system. 

Tile distance between the retroreflector and the airborne laser at time t in: 

d(t, k) = (ZT(t)Z(t) + U(»‘)Tu(k).2ZT(t) (1.2) 

Equation (1.2) docs not completely represent the model for the ALRS ineasurcments since the 
accuracy goals specified for the s>-stcm require that this model be considerably refined. 




(l.l) 


The principal difference between d(t. k) and the ALRS measurement is the refraction incre- 
ment rtt. k); however the model used for ALRS also includes a me.^surcment bias. bj. on the i*^ 
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beam, and a random error, v, which is uncorrelatcd in time or between beams. TIjus the ALRS 
range measurement model is given by 


p(t. k. i) = d(t, k) + nt. k) + bj + vtt. k, i) (1.3) 

where 

t s time tliat measurement is made 
k = retrorellector illuminated at time t 
i s beam used. 

The purpose of the ALRS measurement algorithm is to process the range data (P| , P 2 , . . . , Pn) 
and produce an estimate of the retroretlector imsitions It is clear that this must be 

accompanied, at least implicitly, by estimates of the aircraft positions /.Uy;) at the measurement 
times. t(|. Furthermore, the refraction and the biases must be modellal. and their citects 
compensated. 


The biases are modelled as stationary constants to be estimated, llie refraction increment 
r(t. k) is more difficult to analyze. However, the basic model tor rU. k) is that previously used 
for analysis of satellite laser ranging (5) which was deriveil by (.tardner in ((>) and t7). 

For ALRS purposes, the aircraft height has been included in the model, but since more re- 
cent studies by Gardner (8) have indicated that errors iir r(t, k) do not depend mnni azimuth, 
the c.visting formulation shows refraction compensation errors to depend upon barometric pres- 
sure and the gradient of PTK (i.e., pressure X temperature X coefficient related to lapse time) 
(9. 10). 

Variables which appear in the estimation process are collectively relerred to as the estima- 


tion stale, denoted by X, that is 
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^NXl - 


where 


r^HX \ 


N == 5N^ + Njj + 6 


p^, qj. H respectively, the pressure am! gradient FFK at the retroreflector site 
Njj = the number of iiuiepeiulent laser beams 


Nj^ = llie number of rctroreHector sites 


The inclusion of /(t) will be explained below. 


3.2 The nstimation Process 

An enormous quantity of data can l>e acquired during an ALRS aircraft flight. For a six 
hour flight, the system can operate at a repetition rate of lOpps and receive return pulses on six 
(Njj) beams to register nearly 1,3 million range measurements. This aspect alone rules out any 
“Batch Processing” technique in which a large information matrix is inverted to estimate all 
parameters at once. A sequential estimation process (minimum variance/Kalman filter) has there- 
fore been developed, in which each return, or the returns from a single pulse, are used to 
modify the existing estimate of This approach has been assumed in the definition 

of X^i^x 1 ) id)ove, where a time-varying aircraft pi'»silion appears in the slate, rather than a se- 
quence of independent aircraft positions. In addition, it is assumed explicitly that during a single 
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ALUS fliRlit, tlic other components of the state arc constant; i.c., no crustal motion takes place, 
the beam biases remain the same, ami the error in knowledge ot atmospheric conditions at each 
rctroreflector position remains unchanged. 

At least four distinct concepts are involved in using a sequential algorithm: 

(i) How to change a vector estimate, given a scalar measurement: 

(ii) Modify (i) to include statistical error concepts; 

(iii) How is the process to be initiated; 

(iv) How are the time-varying aircraft positions handled. 

These concepts are treated exhaustively in the literature, (ID among others: thus only a brief 
analysis will be discussed here. 

Tile approach is fundamentally based on linearization. Because a prior estimate exists, an 
estimated or predicted range measurement p(t, k, i) can be computed, based on the estimated 
variables. The actual measurement is then represented in a truncated Taylor’s series as; 

p(t, k. i) = pit. k. i)+-^l ^ (X - X") + V (1.5) 

dX I x = x' 

where a minus or phis superscript red'ers to instant, just before or just after the processing of 
data, respectively. 

This now has the form. 

y = hT(X - X-) + V. (l b) 

One important estimator is the Kalman filter, for which the optimum estimate is expressed 
by 

(X+-X-) = Ky (1-7) 

K being the optimum filter gain. Tlie error, c. in the estimate is 
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c = (X-X-)-Ky = (X-X+) (1.8) 

and the covariance of this error is the expectation of : 


P+ = E(eeT) 


(1.9) 


Since the variances of the error components are given by the diagonal elements of P, 
a minimum variance estimate is obtained by solving for those components of (X - X) which mini- 
mize tr(P) that is: 

tr(P+) = tr(E(ee'f)l 

= tr .^E{X - X-) (X -X‘)T - 2E(X - ^-) yTRT + KE(yyT) 


and 

d(tr P) 
dK 


|-2E(X - X-) yT + 2KE(yyT) j- = 0. 


( 1 . 11 ) 


Then 

K = E(X - X-) yT E(yyT)-t (U 2 ) 

where 

E(X - X-)yT = p-h 

E(yyT) = (hTp-h + Rl 
E(\vT) = R 

Thus 


K = P-hjhTp-h + R)-» 


(1.13) 


The matrix P" is used In (1.13) to obtain the optimum filter gain, which in turn, is used to 
calculate the optimum estimate of the state from (1.7). A new value for this matrix is computed 
from 

P+ = (l-KhT)p- (1.14) 

where I is a unit matrix of proper dimensionality. 


1 .^ 
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Then 

X+ = X- + P-h(hTp-h + R)->hTp-(p - p) (U5) 

The measurement, p, is not exact. There are inaccuracies involved in the pulse timing, model- 
ling errors, etc. All of these error sources are lumped into the uncorrelated noise term (1.5) 
v(t, k, i). On the basis that the primary' source of this error is truncation noise in the clock and 
that bj takes up the constant components, the statistical model for v introduced in (1.5) was 
selected with 

E -^v(t, k, oj- = 0 

E -^v(tn,, k, i) v(tj^, I, = R C *1^ 

where 

gin ^ To m ^ n 
^ \l m = n 

and R is the scalar appearing in (1.12). 

We have shown how a measurement modifies the estimate and its statistics. It is thus ap- 
parent, that it is necessary to provide a starting estimate, and an associated covariance, to initiate 
the process. This has been done by using numbers which can be expected when perfonning the 
actual ALRS process. The apriori positions of the rctrorcflectors will have errors in each com- 
ponent which arc dependent on the care taken in the process of target deployment. 

Tlie aircraft position may have an error standard deviation of 30 to 100 meters because of 
uncertainties in the LORAN/INS position location. Velocity estimates also have a random 
component. 

The estimates of beam biases are zero; preliminary data show that the standard deviation 
about zero is less tlian 1 cm. 
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The estimates of atmospheric parameters have been made by assuming a monitoring station at 
a few retrureflector sites. A standard deviation of 1 mbar is used at the station* 100 mbar elsewhere, 
(i.e., tantamount to no knowledge of the meteorological information at the unmonitored site). 

Extremely large prior covariances could be used, of course, implying heavier weighting of 
the ALRS data. One of the aspects of the data reduction which senes most to increase confi- 
dence in the ALRS approach is that estimation error variances are essentially independent of 
prior variances, provided that a one or more 1 mbar weatlier stations are used. 

In the application of Kalman filtering to systems in which data is taken at discrete time 
points, processing is in two steps: incorporation of the new measurement data and propagation 
of the estimate between measurement points. 

The discussion above has covered how the measurements are used. Between measurements 
all estimates (and true quantities) remain constant except for the aircraft. The aircraft position 
estimate is updated by applying the velocity over the time intersal. For propagation of the co- 
variance, tlie studies reported here have taken tlie conservative stand that each new aircraft posi- 
tion has the same large error variance as does the initial position; thus assuming that the aircraft 
position determined by ALRS cannot be propagated. In actual data reduction, tlic estimated 
aircraft position may be used, lluis improving accuracy; the difficulty of properly evaluating and 
modeling ainrraft disturbances has led us to adopt a more conservative approach for analysis 
purpevses. 

3.3 Coordinates and Constraints 

It is well known (12) that the multilateration problem, of which ALRS is an example, is 
not completely observable; that is, not all of the unknowns in the problem can be determined 
from the ALRS data. The simplest example of this is the fact that the same set of observations 
could be obtained if the complete set of aircraft and retrorefiector positions were translated and 
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rotated as a rigid body. Thus there is a six-fold degeneracy in the problem. To avoid having 
some large variances in the covariance matrix while internal estimates become very accurate - a 
situation leading to numerical problems r an internal coordinate system has been defined for the 
estimation process. In this coordinate system, one rctroreflector is chosen for the origin: one 
retrorcflector defines the x-axis and therefore only its distance from the origin (its x-component) 
is estimated: and finally one retroreflector is taken to define the x-y plane. The variances are 
appropriately modified so that six components arc perfectly known. It is important to note tlut 
the estimates, and variances of the estimates, of baseline length are independent of what coordi* 
nate system is chosen and, except for the numerical problem previously noted, no local coordinate 
system need be chosen. 

These local estimates can be tied back to a larger coordinate system provided independently 
obtained coordinates in the larger system are available for at least three relrorefiectors. 

Baseline data (distance between retrorellectors) is intrinsically local data. When attempting 
to infer subsidence or expansion information it is necessary to define a plane with respect to 
which vertical motion can be measured. The ALRS local coordinate system defines this plane 
with the three “master” retrorcflector locations. - 

4.0 RHSULTS OF SIMULATION STUDIFS 

Simulations have shown that range measurements must be taken at two widely separated 
altitudes in order to strengthen the geometry sufficiently to recover baselines at the centimeter 
level. Thus, in a typical mission, the aircraft approaches the target grid at an altitude of 3.‘>km 
as in Figure 4. After acquiring the first few targets, the instrument shifts to the tracking mode 
for the remainder of the mission. After overtlying the rows of targets at 3.9km, the aircraft 
climbs to its maximum cruise altitude (say 6km) for a second set of passes over the target grid. 
The turning maneuvers between passes can be used to calibrate tlic on board attitude sensors. 
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Figure 4. Typical ALRS Mission Scenario Perspective Showing Ranging from Two Different Altitudes 
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The spacing between targets will depend on several factors including the scientific objectives 
of the mission, the aircraft altitude, and terrain limitations. For the NP3A aircraft, the spacing 
is nominally taken to be 7km. At typical cruise velocities (i.c., 200 knots), and at the most 
favorable aspect angle, laser range data to a given target is taken for approximately 300 seconds 
before the pointing system is commanded to acquire a new target. For a laser repetition rate of 
lOpps, this corresponds to 3000 range measurements per target. Since a given target is common 
to a number of six target sets, and data are collected at two altitudes, over 6500 range measure- 
ments are typically made to each target. 

Figure 5 illustrates the performance of the ALRS evaluated from error analyses. The base- 
line precision vs. baseline distance from an arbitrary origin is shown. The baseline precision de- 
creases with increasing baseline length. For instance, in the absence of atmospheric refraction, 
the baseline precision is 0.65cm for a 20km baseline. The simulation was performed for a 15 
target grid (3 by 5) under the assumption that the single shot laser range measurement uncer- 
tainty was tlcm RMS and the uncorrelatcd biases were on the order of ±lcm. Tlie total number 
of range measurements is 97,959 corresponding to the amount of data collected in approximately 
27 minutes of flight time over the grid assuming no data dropout. 

Tlie baseline precision is degraded slightly in the presence of atmospheric refraction. The 
“with refraction” curve in the Figure was determined under the assumption that surface pressure 
and temperature in the target region could be modelled by quadratic polynomials in the two sur- 
face coordinates and that the coefFicicnts in the polynomials were determined by ground-based 
measurements of pressure and temperature at 15 locations (not coinciding with the target loca 
lions). It was further assumed that the surface measurements of pressure and temperature were 
accurate to ±1 mbar and ±I.4“C respectively. The vertical variation in pressure was assumed to 
be determined by the hydrostatic equation (13). 
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Figure 5. Baseline Free 
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Fif ur. 6 shows ,ho of rrf, action c„o„ and also illns, rates tlntt an catensis. 

network or tneteoroMcal sensors is not rennitej. Attnospheric .neasnrentents made at a sinple 

site ccKlocated with a laser tatjet within tint MRS frid. rednee tlw efreets of 

ntmospherie refraetlon upon AUS baseline nmeision. In the estimation pme-ess desenbed in 
Section 3. if the Infonnation from that sinple ••met sensof is used and tlw atntospintrie param. 
eters for the other ALRS tantets are made pad of the set of estimated parametets. a factor of 
about 7 improsmnent in baseline pmeision is achieved. The also shows that the Inelmion 
of additional meteorological sensors witldn the ALRS tatpet area does not sifninetmtly imptove 
basehne precision. It should be emplutsired that Figure S assumes an extensive meteotolofieal 

sensinp network in the neighborhood of the tarpet repion whereas Fipure 6 utilires meteorolopieal 
sensors collocated at one or more target sites. 

Fipure 7 shows the evolution of basehne pmeision as a function of baseline distance for a 
mpion I4kn. X 1 12km In which 51 laser relroretlectots are deployed. It can be seen ttot the 
baseline precision is depraded at the rate of about I.7en,/I00kn,. This tesult eompams very 
Closely with that obtained for smaller target grid areas. 

In Fipum 8 the evolution of baseline ptecision i, given for a series of randontly deployed 
laser tamets. These tamets are distributed tas shown in the inset) in a potential ALR.S night test 
region in the vicinity of Shenandoah. VA. The targets were located at approximately 17 Hrst 

order survey monuments currently maintained by the U.S. Geological Survey. The alion 

tnd, cates tlut the random tatpet pattern does not signineantly affect the baseline pmeislon tela- 
tive ,0 tlu,, presented in Pigs. 5. 6. and 7. For the simulation, one n.eteorologieal sensor was 

loeatevt in the middle of the grid and the meteotological paran.etets a, other shew wem determined 
in Che cslirnalion process. 
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Figure 6. Airborne Laser Ranging System Baseline Precision as a Function of the Number of Meteorological Stations 
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Figure 7. Baseline Precision vs. Baseline Distance 
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One of the principal applications anticipated for the ALRS surveying system Ls monitoring 
relative motion in the l-arth's upper crust at the sub-centimeter/yr. rate. By developing and 
interpreting the system's ability to detect motions of the Farth's upper crust, a model of the 
strain accumulation aunpatible with ubser\ations of crustal motion and tectonics of a region 
within the ex|>erimental data collection area can be deriwd. 

For example, for a target distribution as shown in I’igunr 7, baseline rate accuracies am be 
obtained from the baseline precision. By assuming that the baselines are measured over a one 
day iH'riod and by repeating the measurements every ninety days, then as sliown in l‘igure 
baseline rate accuracies at the end of a ten year obserx'ation pericn! can be determined for 5 km 
baselines to better than 0.02cm/yr and for 100 km baselines to about 0.1 cm/yr. These baseline 
rates ci>ra*s|H)nd to strain rates of 4 X lO’*' strain/yr and I X UV** strnin/yr. With the baseline 
rate accuracies indicated, the ALHS could provide a capability to obsene the pivcursory geodetic 
motions belie veil to occur before large eartluiuakes. Indeed within a a'gional scale, the Al.KS 
could provide the Unit real pos^ibility of “capturing” a magnitude 7.5 and above eartluiuake. 

CONCLUSIONS 

The Airborne Laser Ranging System is a uniiiue instrument capable of rapidly performing 
dense, large saile geodetic and engineering surveys with subcentimeter accuracies over long base- 
lines. Since all data Is initiated, received, and processed at the aircraft, there is no need for a 
complex data collection network, absolute time information, or skilled personnel in the field. 

It therefore promises to be a highly cost effective device for geophysic's studies, large scale 
surveying, and land management applications. I'sivcially when combined with pliotogrammetric 
instrumentation. Lhis system can detect strain rates to an accuracy of about 5 X 10’^ strain 
per year over a measuring period of 4 years. Such a system can therefore provide the capability 
to observe on a rv'gional scale the precursory crustal motions believed to occur before large 
earthquakes. 
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